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Abstract 


Charge-exchange reactions for scattering to the continuum are 
considered in a high-energy multiple- scattering model. Calculations 
for (p,n) and ( He,H) reactions are made and compared with 
experimental results for 12 C, 16 O, and 27 Al targets. Coherent 
effects are shown to lead to an important role for inelastic 
multiple-scattering terms when light projectiles are considered. 


Introduction 


n high-energy proton collisions, several mechanisms lead to secondary neutron production, 
ow-energy neutrons are produced in evaporation processes. At somewhat higher energies 
secondary neutrons from intranuclear cascades are produced. Neutrons are also produced near 
the incident proton energy from charge-exchange reactions that can occur with or without 
concomitant pion production. These sources of high-energy neutrons may be important for 
risk assessments of space radiation to astronauts. The relative biological effectiveness (RBE) 

° 1S !nnTr g ^ r ) eu J r °“ s . not been studied in great detail; however, one study with neutrons 
near 400 MeV (ref. 1) indicates a biological effectiveness similar to that of 14-MeV neutrons 
with values of RBE between 10 and 30. Because galactic cosmic ray fluxes are dominated by 

high-energy protons, an accurate prediction of secondary neutron sources is then required for 
proper risk assessment. 


In this paper we present a model for neutron production from charge-exchange processes 
a intermediate energies in proton-nucleus reactions. Protons and neutrons may be considered 
as a single particle with an intrinsic degree of freedom called “isospin,” which is analogous 
to the more familiar spin degree of freedom. The charge exchange then corresponds to an 
isospm transition or isospin flip between target and projectile nucleons. Neutron production 
with concurrent pion production (e.g., through isobar formation) will be considered elsewhere. 
From general principles, the free nucleon-nucleon (NN) amplitude for charge exchange (f ex ) is 
written in terms of the proton-proton (pp) and neutron-proton (np) amplitudes as 


Jex — fpp — fnp ( 1 ) 

At high energies (>1 GeV), f pp and f np become approximately equal (ref. 2) and one should 
expect only a small contribution from elastic nucleon events in the nucleus unless nuclear medium 
modifications alter these amplitudes drastically from their free values. We have described a 
model of inclusive (IN) nuclear reactions for high-energy reactions that has been successful in 
escribing momentum distributions for proton and alpha-particle reactions (refs. 3-5). We applv 
this model to the charge-exchange reaction 


p + T — > n + X ^2) 

at intermediate energies where f ex is not negligible. The same model applies to a charge 
exc.ange for neutron projectiles and also to a charge exchange for the mass number A = 3 
nuclei (e.g He + T -> H + X). In the remainder of this paper we present our formalism for 
charge-exchange reactions in inclusive scattering. Physical inputs for our calculations are then 
described, including a description of fits to experimental data for nucleon-nucleon scattering that 
is used to make absolute predictions for nucleon-nucleus reactions. Illustrative calculations of 
the model are then discussed and compared with experiments and, finally, conclusions are made 


Inclusive Scattering Model 

The differential cross section in energy loss and scattering angle for inclusive reactions may 
be written (refs. 3-5) as 

a t 

k 1 f J 2 , ,2 h , ^ q (b-bV(x(b)-X T (b')) 

d£l dEpt ) ^ 


/ip 

d 2 ° \ = J?L_ f d 2 b d %> e^-(b-b') c *(x( b )-x t ( b ')) Wm (b,b',w) (3) 

(2n) 2 J m=l 


After angular integration the energy-loss cross section is 

A t 


= f d 2 b g — 2/mx( b ) jr W m (b,b,w) 
dEpi /in J m= \ 


( 4 ) 


where b and b' are impact parameter vectors, q is the momentum transfer vector fc is the relative 
wave number, X is the eikonal mean-field potential, and the energy loss is u , = E P - E pl , where 
E P and E P , are the initial and final projectile energies. The collision terms in equations (3) 

and (4) are defined as 


1 f m 

W m (b,b',w) = 7-7^2 / II 

(rn!) J j=1 


dk, 


(2 tt)- 


6 ( E f - Ei) < 0 P 0 T |(x (b)) m | Opkj > 


x < kjOp (x 1 (b')) OpOp 


> 


( 5 ) 


where we are considering the ejection of target nucleons with wave number vectors k into 
continuum states, and the ground states of the project, le jmd target are denoted by |0 P > and 
|0 T >, respectively. In equation (5) the eikonal operator (X) is defined as 


c (b) = l Y r dz [ d 3 q e^e- %(i ra e ^ / N N (q) 

(27t) 2 fc N N Q J d-oo j 


( 6 ) 


where a and j label the projectile and target constituents, respectively, tj (b, z) is the 
projectile-target relative separation, r is the internal nuclear coordinate, f m is the nucleon- 
nucleon amplitude, and /cnn is fhe NN relative wave number. 

In references 3-5 we showed that if plane- wave states are assumed for kp the collision terms 
are well approximated by 




Cm £ Cm) 

(m!) J 


711 — 1 


W\ 


( R± ’ S± ’ %)_ 


( 7 ) 


where 


R i = b-b' 


(8) 


and 


S ± = - (b + b r ) 


( 9 ) 


fR is the total binding energy for m knockout nucleons, and the energy-loss parameter is 
o m 

defined as 


2 


U = V 2m N (w-efl m ) (10) 

where wi N is the nucleon mass. In equation (7) the approximation is accurate for forward-peaked 
density matrices, and C x = 1, C 2 = *74, C 3 = */ 105, and C 4 = n 2 / 24lL 

The first collision term is given by 


Wi (b - b>) = r * r f *1 rf v 

(27T) /c^ n 7-oc 7-cx y 




X F (q) F(q') / NN (q) /* N (q') /?i(a,/3,u;) 
where F is the one-body form factor of the projectile, and we have defined 

Q = \ (q + q) 

0 = q - q' 

The target one-particle response function for an uncorrelated wave function is written 
(refs. 3-5) 


( 11 ) 

( 12 ) 
(13) 
as 


#1 («,/3,£i) = 


i 

(27T) 2 


/ 


d 3 a: d 3 y e ia x e ifi y 


Jo(£ix) p (; 


X x\ 
y+ 2' y ~2) 


x0(w-e jBl ) (14) 

where p is the density matrix of the target, j Q is the zeroth-order spherical Bessel function, and 
0 is the unit step function. 

In order to evaluate the inclusive distribution for charge-exchange reactions, we replace 
equation (11) by 


w? (Rj_tSj_,w) = [ dz dz' d 3 a £>(3 e <a ' R e*^' s 

(27r) k z NN J 

x F (a + (3/2) F (a - (3/2) f ex (a + (3/2) f ex (a - (3/2) 

x (a,/3,u;) (15) 

when the projectile loses one unit of charge, for example, (p,n) or ( 3 He, 3 H) reactions. For the 
inverse reactions (n,p) or ( 3 H, 3 He), the factor Z P N T is replaced by N P Z T . In equation (15) 
we are also making the simple approximation that the target response does not distinguish 
between {p,p’) or (p,n) reactions. Also, the energy loss in equation (15) for the (p,n) reaction 
is lu = E P — E n . If we use equation (1) for the two-body inputs, we have 

fex (a + 0/2) f\ x (a - 0/2) = f w (a + 0/2) f ] pp (a - 0/2) 

+ fnp (a 4- 0/2) fl p (a — 0/2) 

-fppioL+0/ 2) fl p (a- 0/2) 

— fpp ( a — 0/2) fnp + 0/2) (16) 
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The contribution from nucleon ejection before or after charge exchange is approximated by 

WT ( ;Ri,Si,ui) = {- 6 W? Wi (Rx,S,^) (17) 

which ignores the noncommutativity of the two-body amplitudes. Higher order processes are 
evaluated similar to equation (17). 


Physical Inputs 

We next describe the physical inputs for our calculations. If we do not consider spin 
dependence, the shell model for a harmonic oscillator basis yields the density matrix 


'(> 


p (y + 2 ,y 2 


x\ _ 1 e -y 2 /R 2 Te - x 2 l AR T 

4) 3/2 

[ ( y 2 x 2 ^ 

x \ a s + 2a p i 2 


flf 4 R l T 


A a d 




(18) 


where Rp is the target matter radii and a s , a p . and a ( i are normalization constants for s-, j>, 
and d-shell nucleons, respectively, given (for A r < 40) by 


4 

Ay 




3 cip — 


i 


At — 4 


12 

Ay 


® a d ~~ i At — 16 


The response function for the Is shell is found as 


( Ay > 4) 
(Ay < 4) 

(Ay £ 4) 

(4 > Ay < 16) 
(Ay > 16) ^ 
(A t < 16) >1 
(Ay > 16) J 


R = e -R 2 T 0 2 /4 e -R 2 T (a 2 H 2 ) sinh 

s (Xy/l! ' 7 


(19) 


( 20 ) 


( 21 ) 


(22) 


for the lp shell as 


R = m N a P^L / [2 - i^ 2 + 2 4 (a 2 + £ 2 )1 Sinh (2i#a*) 

p a0F ([2 v 'J v 


- 4i?^a£ cosh 

-PB,(a 2 + £ 2 ) 
x e T\ ) 





(23) 
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and for the degenerate Id and 25 shells, which we have combined as a single shell because spin 
coordinates are not treated explicitly, 


Rd = 


m N a dRr 

ay/n 


7 - + 2 4 (a 2 + £ 2 ) + 2 4 (<* 2 + £ 2 )' 


+ 8 R^a 2 £ 2 


sinh 


(2 

- 4/3 2 + 4 (« 2 + d)] sinh ( 2 4 aC) 

+ ^4^ 4s ’ n h 

- [ 44 a? + 84«C (a 2 + d)] cosh ( 2 4 «£) 

+ 24/9 2 a^ cosh ( 24 «f) | c -4^/ 4 e "4( a2+ f 2 ) 


(24) 


Values for occupation probabilities calculated from equations (19) (21) and the radius parameter 
Rt are listed in table 1. Separation energies for several nuclei are listed in table 2. (See ref. 6.) 


Table 1. Shell Model Parameters 


Nucleus 

a s 

a p 

® d 

Rt , fm 

4 He 

1.0 

0 

0 

1.33 

6 Li 

.57 

.143 

0 

2.11 

9 Be 

.444 

.185 

0 

1.79 

12 c 

.333 

.222 

0 

1.69 

16 0 

.25 

.25 

0 

1.83 

20 Ne 

.20 

.20 

.033 

2.14 

27 A1 

.148 

.148 

.068 

1.91 

40 Ca 

.1 

.1 

.1 

2.10 


For the 3 He or 3 H projectile we use a form factor F — exp 
( Rp ) of 1.45 fm. The two-body amplitude is parameterized as 


(-$.*) 


( <„\ ct NN (PNN + *) &NN 

/NN (<7) = ^ exp 

4n 


1 


with the projectile radii 


“ ( 2 5 NN^ “ 


l )] 


(25) 


where <tnn is the nucleon-nucleon total cross section, is the ratio of the real part to the 
imaginary part of the forward two-body amplitude, and and C nn are slope parameters. 

Equation (25) is used to fit experimental data (ref. 2) for the elastic pp and up scattering. 
By using equation (1), the charge-exchange data are also fitted with the results plotted as the 
differential cross section versus the Lorentz invariant momentum transfer — t = q 2 , shown in 
figure 1 for several energies. The resulting parameters for equation (25) are listed in table 3. 
The constraint of fitting three sets of data simultaneously does not allow for an accurate fit to 
the charge-exchange data at all values of momentum transfer, especially at lower energies. This 
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Table 2. Binding Energies of Shell Model Orbits 


Nucleus 

Orbit 

£ Sp MeV 

(a) 

4 He 

s 

20.5 

6 Li 

s 

26.0 


P 

5.2 

9 Be 

s 

27.2 


p 

18.2 

12 C 

s 

38.7 


p 

17.5 

1G o 

s 

39.0 


p 

18.0 

27 A 1 

s 

50.0 


p 

22.0 


d 

15.0 

4, »Ca 

s 

51.0 


p 

35.0 


d 

15.0 


a Values for are taken from reference 6. 


Table 3. Nucleon- Nucleon Amplitude Parameters 


T la ,„ MeV 

Ppp 

Pnp 

B pp , fm 2 

B np , fm 2 

C pv , fm 4 

Cfip , fm 4 

460 

-0.9 

0.05 

0.008 

0.1 

0.003 

0.005 

630 

-.54 

-.08 

.23 

.24 

.0065 

.0075 

830 

-.06 

-.35 

.20 

.21 

.0036 

.0036 

1000 

-.05 

-.4 

.22 

.22 

.0036 

.0036 

2200 

-.245 

-.496 

.29 

.282 

.0037 

.0037 


problem, which persists even when more sophisticated parameterizations of the NN amplitudes 
are used, has been attributed to the need for a description of the charge-exchange reaction using 
Quantum Chromodynamics (ref. 7). A more conventional explanation to the behavior of the 
charge-exchange amplitudes than that given in reference 7 is that the amplitude is dominated 
by one pion exchange with large distortion effects from coupling to the elastic channels (ref. 8). 
Because we are not considering interference between the amplitudes in the second-order terms, 
an alternative fit to the charge-exchange data using 


fc.x 


47T 


exp 



(26) 


without regard to equation (1) is also considered; the resulting fits are shown in figure 2 and 
parameters for several energies are listed in table 4. The charge-exchange distributions shown in 
figures 1 and 2 are largely independent of shape. This is in contrast to the elastic pp and np data 
that are nearly isotropic at low energies and show the diffractive forward peak at high energies. 
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Table 4. Charge- Exchange Amplitude Parameters 


Tiab, MeV 

a ex , frrr 

D ex , fm 2 

c,, x , fm 4 

260 

5.9 

0.42 

0.02 

310 

4.27 

.42 

.022 

379 

5.02 

.42 

.02 

460 

3.89 

.36 

.0175 

545 

3.64 

.41 

.017 

630 

3.27 

.46 

.0198 

830 

2.01 

.24 

.0065 

1000 

2.14 

.29 

.007 


Results and Discussion 

We next discuss calculations of charge-exchange differential cross sections for proton and 3 He 
projectiles. In the intermediate energy range (100 to 1000 MeV), three mechanisms occur for 
neutron production in proton-induced reactions. These are (1) neutrons from the intranuclear 
cascade, (2) neutrons from the charge exchange between the incident proton and a bound 
neutron, and (3) neutrons created in the formation and decay of isobars. When comparing 
calculations with experimental results we must keep in mind that only the second mechanism is 
considered in our calculations. Isobar formation becomes important above an incident projectile 
energy of about 500 MeV, and neutrons produced at that level will be peaked at energy losses 
above that corresponding to the pion mass of 139 MeV. Cascade neutrons are those knocked out 
of the target nucleus and will be produced with an energy spectrum peaked at low energies and 
extending out to several hundreds of MeV because of multiple collisions. The neutron spectrum 
from the charge exchange will be peaked at energies close to the incident projectile energy with 
the isospin transitions to bound excited states (not considered here) occurring in the first few 
10’s of MeV. The quasi-elastic knockout, which we will consider, occurs at higher excitation 
energies above the lowest shell separation energies. 

In sketch A we compare calculations for 667-MeV protons on 27 A1 by using the amplitudes 
of equation (25) (solid line) and equation (26) (dashed line). The improved fit of equation (26) 
at smaller momentum transfers (see figs. 1 and 2) leads to a larger peak in the cross-section 
distribution close to the beam energy, whereas at the larger energy losses the results are very 
similar. This comparison demonstrates the importance of having a correct description of nucleon 
charge- exchange amplitudes, especially the slope parameters, in describing the production of 
high-energy neutrons. 

In figure 3, calculations are compared with experimental data for neutron production from 
450-MeV protons incident on an aluminum target of 6.73 g/cm 2 . Figure 3(a) compares the 
experimental data of reference 9 (the hatched area) with the calculations described above (a 
solid line for all order terms and a dashed line for the first-order terms) at a neutron production 
angle of 20°. The histogram shown is the result of the Bertini Monte Carlo model (ref. 9). In 
figure 3(b) a similar comparison is made for a neutron production angle of 30°. The strengths 
of the cross-section distributions are predicted quite well by the calculations; however, the peak 
is shifted to higher energies than the measurements, especially for the comparison at 20°. The 
Monte Carlo results are also in disagreement with the experimental results. One possibility 
for the discrepancy may be due to the thickness of the target (6.73 g/cm 2 ) which presents a 
nonnegligible possibility for multiple collisions of the incident proton or secondary neutron. 

In figure 4 we compare our results with recent experimental data (ref. 10) for 12 C(p, n)X 
reactions at 290 MeV (fig. 4(a)) and 420 MeV (fig. 4(b)) where the cross-section distribution is 
plotted against energy loss (u;). The solid line denotes the calculation described above through 
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10-1 r 


cto/dE, 

mb/MeV 



7, MeV 


Sketch A 


third-order terms and the dashed line denotes a calculation in a relativistic plane- wave model that 
has been normalized to the data as described in reference 10. Agreement between calculations 
and experiment is good except where a small energy loss occurs, which may be due to resonance 
excitation. 

In figure 5, predictions for neutron production from the charge exchange for protons incident 
on 16 0 and 27 A1 are made at several energies. We note that the multiple-shell structure of 
the targets is more apparent at lower incident proton energies because less absorption occurs, 
and this results in a smaller decrease in the contributions from the lower shells in comparison 
with the results at higher energies, i.e., at 1 GeV. At all energies the effects of absorption are 
quite important in reducing the lower shells when compared with the results expected from a 
plane- wave model (ref. 5). 

In figure 6, calculations for the 12 C( 3 He, 3 H)X reaction for an incident 3 He energy of 2 GeV 
are compared with the data of reference 11 for several triton production angles. Agreement at 
the forward angles is quite good on the effects of the quasi-elastic background from the pp and np 
channels, which are seen to be more important than those for incident protons. The agreement 
at larger angles is poor because of the use of an uncorrelated density matrix that is known to 
have insufficient strength at large momentum transfers (ref. 12). 

Concluding Remarks 

A formalism for describing charge-exchange reactions was investigated for light nuclei-induced 
reactions. Comparisons of the calculations with experimental data for target nuclei most 
important for space radiation studies were quite favorable for the model considered. Coherent 
effects were seen to be important for ( 3 He, 3 H) reactions. Improvements in the calculations will 
most likely occur when more realistic charge-exchange amplitudes are considered and by the 
treatment of correlation effects in the nuclear response functions. 


NASA Langley Research Center 
Hampton, VA 23681-0001 
September 17, 1993 


8 



References 

1. Bianchi, Marilena; Baarli, J Sullivan, A. H.; Di Paola, M.; and Quintiliani, M.: RBE Values of 400-MeV 
and 14-MeV Neutrons Using Various Biological Effects. Biological Effects of Neutron Irradiation , International 
Atomic Energy Agency, 1974, pp. 349 357. 

2. Schopper, H., ed.: Elastic and Charge Exchange Scattering of Elementary Particles . Landolt-Bomstein Numerical 
Data and Functional Relationships in Science and Technology , Group I, Volume 7, Springer- Verlag, 1973. 

3. Cucinotta, F. A.; Townsend, L. W.; and Wilson, J. W.: Quasi-Elastic Nuclear Scattering at High Energies. 
NASA TM-4362, 1992. 

4. Cucinotta, Francis A.; Townsend, Lawrence W.: and Wilson, John W.: Multiple-Scattering Effects in Quasi- 
Elastic a - 4 He Scattering. Phys. Review C , 3rd ser., vol. 46, no. 4, Oct. 1992, pp. 1451 1456. 

5. Cucinotta, F. A.; and Dubey, R. R.: Final State Interactions and Inclusive Nuclear Collisions. NASA TP-3353, 
1993. 

6. Barrett, Roger C.; and Jackson, Daphne F.: Nuclear Sizes and Structure. Oxford Univ. Press, 1979. 

7. Fernandez, F.; Valcarce, A.; Gonzalez, P.; and Vento, V.: p(n,p)n and p(p, A ++ )n Charge- Exchange Reactions 
in a Constituent Quark Model. Phys. Lett., vol. B287, no. 1,2,3, 1992, pp. 35-39. 

8. Jain, B. K.; and Santra, A. B.: Study of the p(p, n)A ++ Reaction in Nuclei. Nucl. Phys., vol. A519, 1990, 
pp. 697 720. 

9. Wachter, J. W,; Gibson, W. A.; and Burrus, W. R.: Neutron and Proton Spectra From Targets Bombarded by 
450-MeV Protons. Phys. Review C, vol. 6, no. 5, Nov. 1972, pp. 1496-1508. 

10. Hicks, K. H.; Alford, W. P.; Celler, A.; Henderson, R. S.; Jackson, K. P.; Miller, C. A.; Vetterli, M. C.; Yen, S.; 
Brieva, F.; Horowitz, C. J . ; and Piekarewicz, J.: Comparison of the Quasifree Charge- Exchange Reaction for l2 C 
and 54 Fe. Phys. Review C , vol. 47, no. 1, Jan. 1993, pp. 260 266. 

11. Bergqvist, I.; Brockstedt, A.; Carlen, L.; Ekstrom, L. P.; Jakobsson, B.; Ellegaard, C.; Gaarde, C.; Larsen, J. S.; 
and Goodman, C.: The ( 3 He,Q Reaction at Intermediate Energies— Spin-Isospin Transitions to States in 12 N and 
13 N. Nucl. Phys., vol. A469, no. 4, Aug. 1987, pp. 648-668. 

12. Antonov, A. N.; Hodgson, P. E.; and Petkov, I. Zh.: Nucleon Momentum and Density Distributions in Nuclei. 
Oxford Univ. Press, 1988. 


9 



Experimental data 
□ PP 

O Charge exchange 
Model fits 


— Charge exchange 


-u (GeV/c) 2 



(a) 460 MeV. 


Experimental data 


□ 

PP 


A 

np 


O 

Charge exchange 



Model fits 


— 

PP 


— 

np 

a D 


Charge exchange 


t , (GeV/e) 2 
(b) 630 MeV. 

Figure 1. Comparison of fits using equation (25) with experimental data (ref. 2) for elastic pp 
and up scattering and charge-exchange scattering data. 
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Experimental data 



(c) 830 MeV. 



(d) 1000 MeV. 
Figure 1. Continued. 
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Experimental data 



(e) 2000 MeV. 
Figure L Concluded. 
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(a) 310 MeV. 



(b) 460 MeV. 

Figure 2. Comparison of fits with experimental data (ref. 2) for charge-exchange reactions using 
equation (26). 
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7-iab, MeV 

(a) Neutron production angle of 20°. 



T'lab, MeV 

(b) Neutron production angle of 30°. 

Figure 3. Comparison of calculations with experimental data of reference 9 for 450-MeV protons 
on 27 A1 target of 6.73 g/cm 2 . 
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(a) 290 MeV; neutron production angle of 20.4°. 
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(b) 420 MeV; neutron production angle of 24°. 

Figure 4. Comparison of calculations for 12 C(p, n)X reaction with experimental data of 
reference 10. 
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(c) 1000 MeV. 
Figure 5. Concluded. 
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Figure 6. Comparison of calculations for 12 C( 3 He, 3 H)X reaction at 667 A MeV with experi- 
mental data of reference 11. 
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